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The total r hadronic width can be accurately calculated using analyticity and the operator product expansion. 
The theoretical analysis of this observable is updated to include all available perturbative and non-perturbative 
corrections. The experimental determination of as{M'^) and its actual uncertainties are discussed. 
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1. INTRODUCTION 

The inclusive character of the total r hadronic 
width renders possible an accurate calculation of 
the ratio lil-|7l 



Vr hadrons (7)] 



T\t- 



(1) 



using standard field theoretic methods. If strong 
and electroweak radiative corrections are ignored 
and if the masses of final-state particles are ne- 
glected, the universality of the W coupling to the 
fcrmionic charged currents implies 



Rr ~ N,{\Vud? +\Vus?) 



(2) 



which compares quite well with the experimental 
average Rt — 3.56 ± 0.03. This provides strong 
evidence for the colour degree of freedom Nc- 

The QCD dynamics is able to account quantita- 
tively for the difference between the naive predic- 
tion (H) and the measured value of Rr- Moreover, 
the uncertainties in the theoretical calculation of 
R-T are quite small. The value of Rr can then 
be accurately predicted as a function of as(M^). 
Alternatively, measurements of inclusive r decay 
rates can be used to determine the value of the 
QCD running coupling as{M'^) at the scale of the 
T mass. In fact, r decay is probably the lowest 
energy process from which the running coupling 
constant can be extracted cleanly, without hope- 
less complications from non-perturbative effects. 
The T mass, Mr = 1.7771 ^^ JJ^JfJs GeV @, lies 
fortuitously in a "compromise" region where the 
couphng constant q;s(M^) is large enough that Rr 



is sensitive to its value, yet still small enough that 
the perturbative expansion still converges well. 
Moreover, the non-perturbative contributions to 
the total r-hadronic width are very small. 

It is the inclusive nature of the total semi- 
hadronic decay rate that makes a rigorous the- 
oretical calculation of Rr possible. The only sep- 
arate contributions to Rr that can be calculated 
are those associated with specific quark currents. 
We can calculate the non-strange and strange 
contributions to Rr, and resolve these further 
into vector and axial- vector contributions. Since 
strange decays cannot be resolved experimentally 
into vector and axial- vector contributions, we will 
decompose our predictions for Rr into only three 
categories: 

Rt = Rt.v + Rt.a + Rt.s ■ (3) 

Non-strange semihadronic decays of the r are 
resolved experimentally into vector {Rry) and 
axial-vector {Rr,A) contributions according to 
whether the hadronic final state includes an even 
or odd number of pions. Strange decays {Rr,s) 
are of course identified by the presence of an odd 
number of kaons in the final state. The naive pre- 
dictions for these three ratios are Rr.v — Rt,a — 
{N,/2)\Vud\^ and Rr.s ^ iVc|Ks|', which add up 
to (|). 

2. THEORETICAL FRAMEWORK 

The theoretical analysis of Rr involves the two- 
point correlation functions for the vector V^'^ = 
tpjj^'"ipi and axial-vector A^j — ipjj^''j5ipi colour- 
singlet quark currents ~ u,d,s): 
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t d^xe^'^^0\T{VtAx)VmnO), (4) 
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^ d^xe^'^^{0\T{A^^ix)A'^^{0n0).i5) 



+9^(7" n 



(6) 



The vector (V) and axial-vector (A) correlators 
have the Lorentz decompositions 

(0) 

ij,V/Ay 

where the superscript (J = 0, 1) denotes the an- 
gular momentum in the hadronic rest frame. 

The imaginary parts of the two-point functions 
^iPv/A^l'^^ are proportional to the spectral func- 
tions for hadrons with the corresponding quan- 
tum numbers. The semihadronic decay rate of 
the T can be written as an integral of these spec- 
tral functions over the invariant mass s of the 
final-state hadrons: 



ds 



imn(°)(s) 



(7) 



The appropriate combinations of correlators are 



\Vud\ 



(na.(.) 



n 



(•'^ (s) 



n 



(^) 



(8) 



The contributions coming from the first two terms 
correspond to Rry and Rr^A respectively, while 
i?r,s contains the remaining Cabibbo-suppressed 
contributions. 

Since the hadronic spectral functions are sensi- 
tive to the non-perturbative effects of QCD that 
bind quarks into hadrons, the integrand in Eq. (|^) 
cannot be calculated at present from QCD. Nev- 
ertheless the integral itself can be calculated sys- 
tematically by exploiting the analytic properties 
of the correlators II'^'^^ (s) . They are analytic func- 
tions of s except along the positive real s-axis, 
where their imaginary parts have discontinuities. 
The integral (0) can therefore be expressed as a 
contour integral in the complex s-plane running 
counter-clockwise around the circle \s\ = M^: 



Rr = 67ri 



ds 



2^)n(o-)(.)-2^n(o)(. 



,(9) 



The advantage of expression §) over d) for R^ 
is that it requires the correlators only for complex 
s of order M^, which is significantly larger than 
the scale associated with non-perturbative effects 
in QCD. The short-distance Operator Product 
Expansion (OPE) can therefore be used to orga- 
nize the perturbative and non-perturbative con- 
tributions to the correlators into a systematic ex- 
pansion in powers of 1/s, 

D=2n dhnO=D ^ ' 

where the inner sum is over local gauge-invariant 
scalar operators of dimension _D = 0, 2, 4, . . . The 
possible uncertainties associated with the use of 
the OPE near the time-like axis are absent in 
this case, because the integrand in Eq. in- 
cludes a factor (1 — s/M^)^, which provides a 
double zero at s = M^, effectively suppressing 
the contribution from the region near the branch 
cut. The parameter n in Eq. ( p^ ) is an arbitrary 
factorization scale, which separates long-distance 
non-perturbative effects, which are absorbed into 
the vacuum matrix elements (C(/Lt)), from short- 
distance effects, which belong in the Wilson co- 
efficients C^'^\s,fj,). The D = term (unit oper- 
ator) corresponds to the pure perturbative con- 
tributions, neglecting quark masses. The lead- 
ing quark-mass corrections generate the D — 2 
term. The first dynamical operators involving 
non-perturbative physics appear at Z? = 4. In- 
serting the functions (|l^) into (^) and evaluating 
the contour integral, Rt can be expressed as an 
expansion in powers of with coefficients 

that depend only logarithmically on Mr . 

It is convenient to express the corrections to 
Rr from dimension-D operators in terms of the 
fractional corrections S^^}r,. to the naive contri- 

V^/A 

bution from the current with quantum numbers 
ij, V or ij, A: 



Rt,V/A = -^l^^^dl'^^EW I 1 



EW 



/ . "ud,V/A 
D / 



Rr.s = 3|K.|'W + + (11) 
^'i^^ = ('^li^V' + '^Ij?'!)/^ ^'^^ average of the vec- 



2 



tor and axial- vector corrections. The dimension-0 
contribution is the purely perturbative correction 
neglecting quark masses, which is the same for 
all the components of Rr'- i^^jV/A ~ ^^'^^ ■ 
factors Sew and 5'^^^ contain the known elec- 
troweak corrections. Adding the three terms, the 
total ratio Rt is 

Rr^i{\Vud? + \Vus?) Sew\i + S'ew + '^^"^ 



where the functions 



^ (cos^^c^^+sin^e, 



^ us 



, (12) 



where sin^ Oc = + 



2.1. Perturbative corrections 

In the chiral limit (m„ = md = rus = 0), the 
vector and axial-vector currents are conserved. 
This implies sn'^''^(s) = 0; therefore, only the 
correlator n("+^^(s) contributes to Eq. (^. Ow- 
ing [|lO| to the chiral invariance of massless QCD, 



n: 



(0+ 



...,v (^) = n(°+^\s) ^ n(s)^ {i ^ j) at any fi- 
nite order in a^. Moreover, this result is flavour- 
independent. 

The perturbative QCD contribution to Rr can 
then be extracted from the analogous calculation 
of the ratio R^+e- (*) for e+e^ annihilation, 



Re+e-is) 



cr(e+e — *■ hadrons) 



cr(e+e~ 
= 127rlmnp 



(13) 



where nom(s) is the correlator associated with 
the conserved electromagnetic current Jl^^ = 
Si Qi'4^il^i^i (* = d, s, ...). The result is more 
conveniently expressed in terms of the logarithmic 
derivative of the two-point correlation function of 
the vector (axial) current, 

n=0 ^ ' 

which satisfies a homogeneous Renormalization 
Group Equation. The Kn coefficients are known 



to order 



111-131. For three flavours, one has: 



Kq^K^^ 1; = 1.6398; K^iMS) = 6.3711. 
The perturbative component of Rt is given by 



Rr 



3{l-t-(5(")} 



n=0 



K„A"'\a,), (15) 



ds f asi—s) 



s|=M2 



l-2^ + 277fi -77s ) (16) 



M6 



are contour integrals in the complex plane, which 
only depend on as{AI^). 

The running coupling as{—s) in Eq. (|l^ ) 
can be expanded in powers of as{M^), with co- 
efficients that are polynomials in log (—s/M^). 
The perturbative expansion of S^^"^ in powers 
of ar = as{M^)/Tr then takes [|l|-§ the formQ 
= 'I2n=i(^n + gn)a", whcrc the g„ coeffi- 
cients depend on Kj^^^ and on f3m<n [f^i = —9/2, 
/32 = -8, (iaiWS) = -3863/192, are the coef- 
ficients of the QCD /3-function] : 



19 

K3 - —K2P1 



19 

24 



265 - 247r^ 



288 



Plj al + O(a^) 

= ar + 5.2023 al + 26.366 al + 0«). (17) 

One observes ||] that the gn contributions are 
larger than the direct contributions ((72 = 
3.563, (73 = 19.99). For instance, the bold-guess 
value K/^ ^ K^{K^/ K2) ~ 25 is to be com- 
pared with 34 = 78. These large "running" con- 
tributions are responsible for the sizeable renor- 
malization scale dependence found in ref. [p^ . 
The reason of such uncomfortably large contri- 
butions 1^ stems from the complex integration 
along the circle s — exp {i<j>) (0e[0, 27r]) in 
Eq. (|l^), which generates the gn terms. When 
the running coupling as(— s) is expanded in pow- 
ers of as(M^), one gets imaginary logarithms, 
log (— s/M^) = i((^ — tt), which are large in some 
parts of the integration range. The radius of con- 
vergence of this expansion is actually quite small. 
A numerical analysis of the series |5) shows that, 
at the three-loop level, an upper estimate for the 
convergence radius is ar,conv < 0.11. 

^ In ref. the perturbative contribution to Rt was 

parametrized in terms of the coefficients F„ , appearing in 
the expansion of the spectral function Imn(s) in powers 
of as{s)/TT. Both parametrizations are related by trivial 
factors: K2 = F3; K3 = F4 + {n'^ /Sf / 12) . 
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Note|^ however, that there is no deep rea- 
son to stop the yl^")(Q;i) integral expansions at 
0{a^). One can calculate the A^"-\as) expan- 
sion to all orders in as, apart from the unknown 
f3n>3 contributions, which are likely to be negli- 
gible. Even for Ur larger than the radius of con- 
vergence ar.conv, the integrals A^"-\as) are well- 
defined functions that can be numerically com- 
puted, by using in Eq. (|l^) the exact solution for 
as{s) obtained from the renormalization-group 
function equation. Thus a more appropriate ap- 



proach 
in Eq. 



is to use a Kn expansion of as 



(15|), and to fully keep the known three- 
loop-level calculation of the functions A("^(q;s). 
The perturbative uncertainties are then reduced 
to the corrections coming from the unknown Pn>3 
and Kn>3 contributions, since the g„ contribu- 
tions are properly rcsummcd to all orders. To ap- 
preciate the size of the effect. Table |l| gives the ex- 
act resuhs [|l for A^'^^as) (n = 1,2,3) obtained 
at the one-, two- and three-loop approximations 
(i.e. /3„>i = 0, /3„>2 = 0, and /3„>3 = 0, respec- 
tively), together with the final value of (5'"', for 
Ct- = 0.1. For comparison, the numbers coming 
from the truncated expressions at order are 
also given. Although the difference between the 
exact and truncated results represents a tiny 0.6% 
effect on Rr, it produces a sizeable 4% shift on the 
value of S^^\ The J'-*'-' shift, which reflects into a 
corresponding shift in the experimental as(-M^) 
determination, depends strongly on the value of 
the coupling constant; for = 0.14 the (J^"-* shift 
reaches the —20% level. 

Notice that the difference between using the 
one- or two-loop approximation to the /3-function 
is already quite small (1.4% effect on J^°^), while 
the change induced by the three-loop corrections 
is completely negligible (0.1%). Therefore (unless 
the /9-function has some unexpected pathological 
behaviour at higher orders), the error induced by 
the truncation of the /3-function at third order 
should be smaller than 0.1% and therefore can be 
safely neglected. 

The only relevant perturbative uncertainties 
come from the unknown higher-order coefficients 
Kn>3- To obtain an estimate of the error induced 



Table 1 

Exact results for A^^^Us) {n — 1,2,3) obtained 
at the one- (/?„>i — 0), two- (/3„>2 = 0) and 
three-loop (/3„>3 — 0) approximations, together 
with the final value of S^^\ for Ut — 0.1. For 
comparison, the numbers coming from the trun- 
cated expressions at order a'^ are also given 



Loops 


^(1) 








1 


0.13247 


0.01570 


0.00170 


0.1690 


2 


0.13523 


0.01575 


0.00163 


0.1714 


3 


0.13540 


0.01565 


0.00160 


0.1712 




0.14394 


0.01713 


0.00100 


0.1784 



Table 2 

^(0) fQj. different values of = asiM^ 



as 


^(0) 


as 


^(0) 


0.24 


0.118± 0.003 


0.34 


0.191 ±0.009 


0.26 


0.132 ±0.004 


0.36 


0.205 ±0.010 


0.28 


0.146 ±0.005 


0.38 


0.220 ±0.012 


0.30 


0.161 ±0.006 


0.40 


0.234 ±0.013 


0.32 


0.176 ±0.008 


0.42 


0.248 ±0.013 



A similar suggestion was made in ref. ]l5|. 



on 5(0), we win take A(5(")) - ±K4A('^\as). 
The naive guess K4 ~ [K^/ K2)Kj, « 25 gives, 
for ar = 0.1, a small A((5(°)) = ±0.004 effect. 
The sensitivity on the choice of renormalization 
scale and renormalization scheme has been stud- 
ied in ref. Q , where it has been shown to be very 
small. 

The perturbative contribution (S'"' , obtained 
through Eqs. (|l5|) and (|l|), is given in Table H 
for different values of the strong coupling con- 
stant q;s(M^). In order to be conservative, and 
to account for all possible sources of perturbative 
uncertainties, we have used Q ~ 50 for esti- 
mating A((5('')). 

2.2. Leading quark- mass corrections 

The contributions (5|^' to the ratio Rt are 

simply the leading quark-mass corrections to the 
perturbative QCD result of the previous section. 
These contributions are known [^,^ to order . 
Quark-mass corrections are certainly tiny for the 



4 



up and down quarks ((5^^ ~ —0.08%), but the 
correction from the strange-quark mass is impor- 
tant for strange decays [p[p^: (niu = rud = 0) 

-2 



5(2) 



+ + 11.03 



(18) 



where rus = ms{M^) is the running mass of the 
strange quark evaluated at the scale Mr- For 
as{M^) = 0.33, 6^uJ «^ -19%; nevertheless, be- 
cause of the sin^ dc suppression, the effect on the 
total ratio Rr is only -(0.9 ± 0.2)%. 

2.3. Non-perturbative contributions 

Since the t mass is a quite low energy scale, 
we should worry about possible non-perturbative 
effects. In the framework of the OPE, the long- 
distance dynamics is absorbed into the vacuum 
matrix elements {0(fi)), which are (at present) 
quantities to be fixed phenomcnologically. If the 
logarithmic dependence of the Wilson coefficients 
C^'^\s,fj,) on s is neglected (this is an effect of 
order al), the contour integrals can be evaluated 
trivially using Cauchy's residue theorem, and are 
non-zero only for D = 2,4,6 and 8. The correc- 
tions simplify even further if we also take the chi- 
ral limit {rriu = md = rris = 0). The dimension- 
2 corrections then vanish because there are no 
operators of dimension 2. In the chiral limit, 
sn(°)(s) = 0; thus only the n(°+i)(s) term in 
Eq. (^ contributes to Rr- The form of the kine- 
matical factor multiplying n(''+^)(s) in Eq. (^ is 
such that, when the s-dependence of the Wilson 
coefficients is ignored, only the D = 6 and D = 8 
contributions survive the integration. The power 
corrections to R^ then reduce to 



-247r^ 



M6 



EC 



(0+1) 

ijy/A 



(O) 



(D=8 



Ml 



(19) 



(20) 



and ^IJ";)/^ ~ for 2n^6,i 

When the logarithmic dependence of the Wil- 
son coefficients on s is taken into account, oper- 
ators of dimensions other than 6 and 8 do con- 
tribute, but they are suppressed by two pow- 
ers of as{M^). The largest power corrections to 



Rr then come from dimension-6 operators, which 
have no such suppression. Their size has been es- 
timated in ref. Q, using published phenomeno- 
logical fits to different sets of data: 



-(0.7 ±0.4)%. (21) 



These power corrections are numerically very 
small, which is due to the fact that they fall off 
like the sixth power of l/AIr- Moreover, there 
is a large cancellation between the vector and 
axial-vector contributions to the total hadronic 
width (the operator with the largest Wilson coef- 
ficient contributes with opposite signs to the vec- 
tor and axial-vector correlators, due to the 75 
flip). Thus, the non-perturbative corrections to 
Rr are smaller than the corresponding contribu- 



tions to Rr 



y/A- 



A more detailed study of non- 



perturbative corrections, including the very small 
D = A contributions proportional to quark masses 
or to as{M^)^ , can be found in ref. 0. 

2.4. Electroweak corrections 

The electroweak corrections to the ratio Rr are 
quite sizeable, because the corrections to the nu- 
merator include logarithms of Mz /Mr , which are 
not present in the corrections to the denomina- 
tor. These logarithms represent a short-distance 
correction to the low-energy effective four-fermion 
coupling of the r to lyrdu or VrSU, and, therefore, 
should be absorbed into an overall multiplicative 
factor Sew- Using the renormalization group to 
sum up higher-order a"\og^{Mz/Mr) contribu- 
tions. Sew becomes |l3 



Sew — 



«("^b) 
a(M2) 

1.0194. 



a{ml) 



a(Ml) 



(22) 



The residual non-logarithmic electroweak correc- 
tion is quite small |18|] : 

= (5/12) a(Af2)/7r ~ 0.0010. (23) 

3. DETERMINATION OF a^iM^) 

The final QCD predictions ^-01 for Rry, Rr.A, 
Rr,s and Rr are given in Tables^ and ^, as func- 
tions of the coupling constant as{M^)- 
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Table 3 

QCD predictions for the different components of 
the T hadronic width as function of cTg = as{M^) 

c^s Rt,v Rt,a Rt,S 



dent determinations of Rr are 

i?f = 3.56 ±0.04, (26) 
= 3.55 ±0.06. (27) 

Using the predictions in Table ^, the average of 
the two experimental determinations of Rr, 

i?^ = 3.56 ±0.03, (28) 

corresponds to 

as (M^) 0.33 ± 0.03. (29) 

Once the running coupling constant Usi^J.'^) is 
determined at the scale Mi-, it can be evolved to 
higher energies using the renormalization group. 
The error bar on as(^^) must also be evolved 
using the renormalization group. Its size scales 
roughly as as(^^)^, and it therefore shrinks as /x 
increases. Thus a modest precision in the deter- 
mination of as at low energies results in a very 
high precision in the coupling constant at high 
energies. After evolution up to the scale Mz, the 
strong coupling constant in Eq. ( p9| ) decreases to 

a.(Af|) = 0.120 1°:™^ (30) 

in excellent agreement with the present LEP av- 
erage (without Rr) a,(M|) = 0.123± 0.006 
and with a similar error bar. The comparison 
of these two determinations of q;s(/z^) in two ex- 
treme energy regimes, Mt and AIz, provides a 
beautiful test of the predicted running of the 
QCD coupling constant. 

3.1. Semi-inclusive r-decay widhts 

With a^iM^) fixed to the value in Eq. (H), 
Table ^ gives definite predictions for the semi- 
inclusive T-decay widths: Rt,v = 1-76 ± 0.06, 
i?^^ = 1.66±0.06, Rr,s = 0.145±0.006. The ex- 
perimental results on exclusive hadronic r-decay 
modes provide then a consistency check of the re- 
liability of the theoretical analysis. 

The assignment of a given measurement to one 
of the three categories Rt,v, Rt,a, Rt,s, is not 
completely straightforward. One needs to have a 
clean tt/K identification, and to know the exact 
number of neutral particles, to separate the vector 
and axial-vector contributions. Fortunately, the 
improved quality of the recent data has allowed 



0.24 


1 


66 ±0 


02 


1 


56 ±0 


03 





145 ±0 


005 


0.26 


1 


68 ±0 


02 


1 


58 ±0 


03 





145 ±0 


005 


0.28 


1 


70 ±0 


02 


1 


61 ±0 


03 





145 ±0 


005 


0.30 


1 


72 ±0 


02 


1 


63 ±0 


03 





145 ±0 


006 


0.32 


1 


75 ±0 


02 


1 


65 ±0 


03 





145 ±0 


006 


0.34 


1 


77 ±0 


02 


1 


67 ±0 


03 





145 ±0 


006 


0.36 


1 


79 ±0 


02 


1 


69 ±0 


03 





144 ±0 


006 


0.38 


1 


81 ±0 


03 


1 


71 ±0 


03 





144 ±0 


007 


0.40 


1 


83 ±0 


03 


1 


73 ±0 


03 





143 ±0 


007 


0.42 


1 


85 ±0 


03 


1 


75 ±0 


04 





143 ±0 


007 



Table 4 

Rr for different values of = as(M^) 



as 


Rr 


as 


Rr 


0.24 


3.37 ±0.02 


0.34 


3.58 ±0.03 


0.26 


3.41 ±0.02 


0.36 


3.63 ±0.03 


0.28 


3.45 ±0.02 


0.38 


3.67 ±0.04 


0.30 


3.50 ±0.02 


0.40 


3.71 ±0.04 


0.32 


3.54 ±0.03 


0.42 


3.75 ±0.04 



The experimental value for Rr is actually deter- 
mined by measuring the leptonic branching frac- 
tions: 



R' 



1 — Be — 

b: 



(24) 

Vri^ i^fXl)]/^T and Tr is the 



where Bg = T[t~ 
total decay rate. 

An independent determination of Rr can be 
obtained by measuring the lifetime. Because the 
decays — > Vri~vi{j) are purely electroweak 
processes, their rates Tr^e = Bi/xr can be calcu- 
lated theoretically with great accuracy. The only 
unknown in Eq. ( p^ ) is therefore the total decay 
rate H: 



„r (1.632 ± 0.002) X 10- 

R^ = 



1.97257. (25) 



The present p| results for these two indepen- 
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to perform an explicit identification of multiple 
7r°'s |2l|,^ and the first systematic analysis 
of exclusive Cabibbo-suppressed decays. 

The Particle Data Group ^ has made a con- 
strained fit to the data, using 12 basis modes 
whose branching fractions sum exactly to unity: 
e~Vf.VT^ 7r~z/^, K~Vr, if *(892)~i^T-, 

TT^ir'^i/r, h^2T:^VT, h^iir^VT, h^Air^Vr, ^h^h^i/r, 
2h-h+{> lh°)iyr and 3ft.-2/i+(> Oh°)iyr, where 
and denote charged and neutral hadrons, 
respectively. Assigning the modes with an un- 
known number of neutral pions to the category 
corresponding to the minimum possible multi- 
plicity [i.e. the decay into 2/i^/i+(> to 
Rry and 3ft.~2/i+(> 0/^°)^'^- to Rt,a], this fit pro- 
vides a first approximation to the actual semi- 
inclusive decay widths. The resulting values can 
then be slightly corrected using the more exclu- 
sive measurements of the modes 3h^ 2h^ tt'^ 1/^- , 

2h-h+2TT°Ur, m^TT^Vr, K*{m2fTT-{> Oh°)Ur, 

K-K+TT-{> Oh°)ur and fs:*(892)-(> Oh^)vr- 
Following this procedure and taking into account 
the most recent data |2|-|24), I geiQ: 

Rry = 1.76 ±0.04, 

= 1.66 ±0.05, (31) 
R^s = 0.14 ±0.02, 

which add to Rr = 3.56 ± 0.07. The agreement 
with the theoretical predictions is excellent. 

The Rr,s predictions are very sensitive to the 
power corrections. As shown in Table there 
is practically no dependence on the value of the 
strong coupling in this case. In fact, the final pre- 
dictions turn out to be very close to the nai've ex- 
pectation Rr^s ~ A^dKisP — 0.147, because there 
is a strong cancellation between the perturbative 
contribution 5^^'^ and the strange-quark-mass cor- 
rection (J^^-* . The success of the theoretical pre- 
diction can then be taken in this case as a test of 
the D = 2 contribution. With more precise data, 
it could be possible to use i?T,s to check the value 
of the strange-quark mass. 

3.2. e+e test 

In the vector channel, one can also use the in- 
formation obtained, through an isospin rotation, 

* The Hi r en t sum of the measured exclusive modes into 
kaons gives Rt-,s = 0.16 ± 0.02. 



from the isovector part of the e^e~ annihilation 
cross-section into hadrons. The exclusive r-decay 
width into Cabibbo-allowed modes with = 1~ 
can be expressed as an integral over the corre- 
sponding e+e~ cross-section p5| , 

3C0S^ 0c ^ f^'^^ , ,..2 ^2 

X (Af 2 + 2s)sY, e- ^yo (s) . (32) 

The analysis of the separate exclusive vector 
modes [^,^ shows a good agreement between 
the actual r-decay measurements and the num- 
bers obtained from e~^e~ data. However, the 
r-decay data is already more accurate than the 
e+e~ results; therefore, this exercise does not im- 
prove the Rt,v measurement in Eq. (|3^). Never- 
theless, the e+e" data offers us the opportunity to 
make an additional test of the QCD predictions, 
by varying the value of M = Mr in Eq. ( ^ . 

The theoretical predictions for Rr,v as a func- 
tion of M can be trivially obtained from the 
formulae given in refs. For lower values 

of M it is important to use the resummed per- 
turbative expansion of ref. ||], since the bigger 
values of as(Af^) in that region imply that the 
non- resummed expansion is non-convergent. Fig- 
ure H^] compares the theoretical predictions 
with the results obtained from e+e" data (the 
points with vertical error bars). Above 2 GeV the 
data is rather conflicting and incomplete, and one 
needs to rely on extrapolations of the fits done at 
lower energies; so this region has not been plot- 
ted]]. The shaded area between the two dashed 
curves corresponds to the theoretical prediction 
for as{M^) = 0.33. The big allowed region at 
low values of M is due to the uncertainty in the 
leading non-perturbative correction, taken from 
the estimate in Eq. (|2l|), i.e. 

S'^^M) = (0.024 ± 0.013) X [Mr/Mf . (33) 

At the r-mass scale, Sy\Mr) = Sy^ is a very 
tiny correction; however, since it scales as the 

^ Although one finds also quite good agreement between 
theory and experiment above 2 GeV, the big error bars 
of the "experimental" points make the comparison quite 
meaningless. 
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Figure 1. Rt,v as a function of M = Mr- 



sixth power of M, this non-perturbative contri- 
bution (and its associated uncertainty) increases 
very rapidly as M decreases. 

Allowing the value of as(M^) to change in the 
range , one gets the larger region between the 
two dot-dashed curves in Figure I. For M ~ Mr, 
the theoretical uncertainty is dominated by the 
uncertainty in the input value of as(M^); the er- 
ror on S'y' becomes dominant for i\/ < 1.6 GeV, 
and overwhelms the result for M < 1.2 GeV. 

One can notice that there is a good agree- 
ment between the QCD predictions and the e+e~ 
data points for M > 1.2 GeV. This confirms the 
role of the threshold factor (1 — s/M'^Y, which 
minimizes the theoretical uncertainties near the 
physical cut, and further supports the theoretical 
framework used to analyze the r hadronic width. 

The departure of the theoretical prediction 
from the data points below 1.2 GeV signals the 
important role of higher-order power corrections 
in this region. The subleading dimension-8 cor- 
rection has been neglected before because, at 
the r-mass scale, its contribution is expected to 
be smaller than the uncertainty on 5y\ How- 
ever, when going to smaller values of Af, the 
5y {M) contribution increases much faster than 
the dimension-6 one and at some point would 
even become dominant, indicating a breakdown of 



the expansion in powers of 1/M . We can use the 
lower-mass data points to make an estimate of the 
size of this contribution. Taking as(M^) = 0.33 
and 5y^ — 0.024, a quite reasonable fit is obtained 
for (5y ' = —0.0095. This is shown by the continu- 
ous curve in Figure 0. Although the S^iM) cor- 
rection is tiny at M = Mr, its effect changes com- 
pletely the predicted behaviour below 1.2 GeV. 
Note, however, that for this value of 5y^ one has 
5'^\m) / 5'^\m) = -1.25 at M = 1 GeV, which 
puts some doubts on the applicability of the in- 
verse power expansion at such a low scale. If one 
takes only into account the region above 1.2 GeV, 
the size of the experimental error bars does not 
allow us to make a clear statement about the size 
of 6y {5y = is compatible with the data), al- 
though smaller values of 6y seem to be preferred. 

3.3. Hadronic invariant-mass distribution 

The leading non-perturbative contributions to 
Rr are suppressed by a factor of 1/M^ and, there- 
fore, are very small. Nevertheless they introduce 
a small uncertainty in the predictions, since their 
actual evaluation involves a mixture of experi- 
mental measurements and theoretical consider- 
ations, which are model-dependent to some ex- 
tent. It would be better to directly measure this 
contribution from the r-decay data themselves. 
This information can be extracted j2|] from the 
invariant-mass distribution of the final hadrons in 
r decay. 

Although the distributions themselves cannot 
be predicted at present, certain weighted integrals 
of the hadronic spectral functions can be calcu- 
lated in the same way as Rr- The analyticity 
properties of H^^^^^^^^ imply |2|,||]: 

dsw{s)iv^n'^;l^y^^ = 

'^'^^'^^^v/A^ (34) 

with W{s) an arbitrary weight function without 
singularities in the region \s\ < sq. Generally 
speaking, the accuracy of the theoretical predic- 
tions can be much worse than the one of Rr, be- 
cause non-perturbative effects are not necessar- 
ily suppressed. In fact, choosing an appropriate 
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weight function, non-perturbative effects can even 
be made to dominate the final result. But this is 
precisely what makes these integrals interesting: 
they can be used to measure the parameters char- 
acterizing the non-perturbative dynamics. 

To perform an experimental analysis, it is con- 
venient to use moments of the directly measured 
invariant-mass distribution {k,l > 0) 



ds 1 



So 



\ ] (35) 



M. 



ds 



The factor (1 — s/sq)*^ supplements (1 — s/M^)"^ 
for So 7^ A/^, in order to squeeze the integrand at 
the crossing of the positive real-axis and, there- 
fore, improves the reliability of the OPE analysis; 
moreover, for sq = it reduces the contribu- 
tion from the tail of the distribution, which is 
badly defined experimentally. A combined fit of 
different i?^'(so) moments results in experimental 
values for as{M^) and for the coefficients of the 
inverse-power corrections in the OPE. 

The first experimental study has been done 
by ALEPH @, using the moments (0,0), (1,0), 
(1,1), (1,2) and (1,3). R°° = Rr uses the over- 
all normalization of the hadronic distribution, 
while the ratios D^'- = R'f / Rt are based on the 
shape of the s distribution and are more depen- 
dent on non-perturbative effects 1 29 1 . The result- 
ing strong coupling constant measurement [ ^0| , 
as{M'^) = 0.330±0.046, is in excellent agreement 
with the Rr determination in Eq. (|29|). More- 
over, the ALEPH analysis has determined the to- 
tal non-perturbative contribution to Rt to be 



r^NP = (0.3 ± 0.5)%, 



(36) 



confirming the predicted ^ suppression of non- 
perturbative corrections. 

More recently, ALEPH Q and CLEO H 
have reported preliminary results of two inde- 
pendent analyses, with a much larger data sam- 
ple. Although their measured s-distributions are 
in excellent agreement |3ll , the fitted Us values, 
shown in Table ||, disagree by about 2 standard 
deviations. The origin of the discrepancy stems 
from their different normalizations. While CLEO 
uses the world-averaged Bi values quoted by the 
Particle Data Group B , ALEPH has smaller (and 



Table 5 

Fitted values of q;s(M^) and the non-perturbative 
corrections. {^GG) is given in units of GeV^ 





ALEPH 


CLEO 


a.(M2) 


0.355 ±0.021 


0.309 ±0.024 


i^GG) 


0.006 ±0.012 


0.024 ±0.010 


5(6) 


0.006 ±0.006 


0.014 ±0.004 


_5(8) 


0.0017 ±0.0014 


0.0028 ±0.0006 


xVd.o.f. 


0.002 


2.1 



precise) leptonic branching ratios and therefore a 
larger r hadronic width, Rr = 3.645± 0.024. The 
average of the ALEPH and CLEO asiM^) deter- 
minations, 



as{M^) = 0.335 ±0.016, 



(37) 



nicely reproduces the result (P9[), obtained from 
Rr using the theoretical estimates of the small 
non-perturbative corrections. Nevertheless, the 
error in (|37| ) should probably be increased by a 
factor of 2, to account for the present 2a discrep- 
ancy between the individual measurements. 

In addition to confirm that the total non- 
perturbative correction to Rr is indeed small, 
the moment analysis allows to make a determi- 
nation of the D — 4 gluon condensate (which 
does not play any role in Rr), in good agree- 
ment with the previous phenomenological esti- 
mate i^GG) = 0.02 ±0.01 GeV^ |,|. 

4. THEORETICAL UNCERTAINTIES 

The accuracy of the experimental data is obvi- 
ously going to improve in the near future. Thus, 
the final error in the determination of as (M^) will 
be limited by the accuracy of the theoretical cal- 
culation. It is then important to carefully study 
the possible sources of theoretical uncertainties 
and try to pin down their numerical effects. 

4.1. Higher-order perturbative corrections 

The precision of the perturbative calculation 
of Rr has been extensively discussed in ref. 
Once the important higher-order logarithmic ef- 
fects have been properly identified and resummed 
to all orders, the size of the remaining renormal- 
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ization scheme (and scale) dependence is quite 
small. The only relevant perturbative uncertain- 
ties come from the unknown higher-order coeffi- 
cients Kn>3- We have estimated their numerical 
effect by assuming an algebraic growth of these 
coefficients and, moreover, in order to be conser- 
vative, we have further increased this estimate by 
a factor of two @, i.e. - 2{Ks/K2)Ks ~ 50: 

A((5(*')) =±i^4A('')(a,) =±0.008. (38) 

Clearly, this is an arbitrary prescription, but it is 
probably the most conservative one in the absence 
of a direct calculation of i4'4. Note, that this value 
of K^^ corresponds already to a non-convergent 
perturbative series: 

K^A^^\as)=QmQr^ KiA^^\as). (39) 

A recent estimate of Ki^^ based on the application 
of scheme-invariant methods such as the principle 
of minimal sensitivity |33|] or the effective charges 
approach [Q, finds ||35| = 27.46. A similar 
value Kj^ = 29 ± 4 ± 2 has been obtained from a 
direct fit to the experimental data The two 
estimates are indeed smaller than our conserva- 
tive choice Ki — 50. 

4.2. Renormalons 

Possible uncertainties related to the asymptotic 
(at best) nature of the QCD perturbative expan- 
sions at large orders have been considered re- 
cently [p7|"^ . The leading large-order contribu- 
tions come from the so-called renormalons, which 
are associated with known singularities p8|] in the 
Borel transform 

1 ft" 
B[D]^b{h)^—^Y.K^^^^-^ (40) 

n=0 

of D{a) = X^n ^nfl"/(47''^), where a = Us/'k- In 
terms of Feynman graphs, the renormalons corre- 
spond to inserting a chain of quark-loop bubbles 
into a gluonic line. Diagrams of this kind gen- 
erate a factorial growth of the perturbative D{a) 
series (i.e. ^ n\) at large orders. Owing to 
the explicit l/n\ factor in the coefficients of (^), 
D{h) is a much better behaved series. Thus, one 
could hope to define D{a) through its Borel sum 

/•OO 

D{a)-D{Q)^l dbD{b) cxp{-b/a). (41) 

^0 



Unfortunately, the function D(b) has singularities 
at 6 = mbo = 2m/(-/3i) (m = -1, ±2, ±3, ...) 
and therefore the integral ( p4| ) is not defined. 

The singularities at m — -1-2, -1-3, ... (infrared 
renormalons) are generated by the low-momenta 
behaviour of these higher-order diagrams, which 
produces contributions of the form D{a) ~ 
X;„?^! (a/"^^o)"; i.e. b{b) ^ -mbo/{b - mbo) 
for b < mbo- The pole at 6 = mbo gives rise to 
an ambiguity when one tries to reconstruct D(a) 
from b{b): [1/a w log (^/i/A)] 

SD{a) - exp i~mbo/a) - (A^s)"' . (42) 

These infrared 1 / s™ contributions are reabsorbed 
p8|-^0|] into the non-perturbative terms of the 
OPE. 

The absence of a singularity at m = 1 is in 
fact related to the lack of any physical gauge- 
invariant local operator of dimension 2. This 
was challenged in ref. |^l|, where the relation 
6[Imn](6) - [sininb/bo)/ i7rb/bo)]bib) was es- 
tablished, implying that either D{b) has a pole 
corresponding to to = 1 or ;B[Imn](6) has a zero 
at that point. The first possibility (an infrared 
renormalon for to = 1) would imply a 1/s ambi- 
guity which could not be reabsorbed by the OPE. 
This puzzle has been clarified by the exact calcu- 
lation of i)(6) (both in QED and in QCD ||]) 
in the limit of a large number of quark flavours, 
which shows that there is no infrared renormalon 
corresponding to to = 1, and y8[Imn](6) has in 
fact a zero. This confirms the relation of infrared 
renormalons to gauge-invariant operators |3^Q . 

The ultraviolet renormalons (to = —1,...) are 
generated by the high-momenta behaviour and, 
because of asymptotic freedom, give rise to a 
Borel-summable series (they are in the negative 
real axis, outside the integration region), which 
makes them harmless. They don't put any real 
limit |3^ ] to the applicability of perturbation the- 
ory. Nevertheless, if the Borel summation is not 
performed, the ultraviolet renormalons induce an 
intrinsic uncertainty in the truncated perturba- 
tive series. 

The factorial growth of the perturbative se- 
ries is indeed dominated by the contribution of 
the leading ultraviolet renormalon (to = — 1), 
Kn ^ n!/(— 6o)"j which gives rise to an asymp- 
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totic series: 



l^ial > l^sal > ... > \KN-ia 



Af-l| 



(43) 

The successive terms decrease until N ~ bo /a, 
where the minimum value is attained, and the 
series explodes afterwards. The alternating sign 
of the Kn coefficients guarantees that the series is 
Borel-summable. However, if one only considers 
the truncated series at a given finite order, the 
accuracy is obviously limited by the size of the 
minimum term MM: 



~ exp i-bo/a) - A^/s . (44) 

Thus, the m — ~1 ultraviolet renormalon 
seems to induce an uncertainty in the perturba- 
tive series, which scales like 1/s. This has been 
used 1^ , ^ to advocate for an additional (small) 
uncertainty in i?r proportional to such a 

conclusion is incorrect, because: 



1. As it stands, the estimate (^4|) depends on 
the chosen renormalization scheme (the def- 
inition of A). A more careful analysis 
shows that the ambiguity actually scales 
as Ay/ajJ^A^s/fi'^, where the scheme de- 
pendence has been absorbed into the factor 
A. The result (|4^ ) only follows if one takes 

= s. Keeping an arbitrary renormal- 
ization scale ^, the A^s//i^ term does not 
contribute to the close-contour integral (||) 
defining R^- (Cauchy theorem). 

2. One could still argue |^ that the improved 
perturbative series, which resums higher- 
order logarithms through the A^") (a^) func- 
tions in Eq. (^6|) , takes jj? = s and thus the 
1/Af^ ambiguity comes back. However, it 
was shown in ref. ||] that the scale depen- 
dence of the resummed perturbative result 
is tiny; taking /i — ^y/s, the difference be- 
tween ^ = 1 and ^ = 3 is very small. A 
A^/(f^Af^) term would instead change by 
a factor of 9, making such a contribution 
irrelevant for ^ ^ 3. 

3. Eq. ( ^4| ) is just an estimate of the ultimate 
perturbative error, defined as the size of the 



minimum term. At the Mt scale, a « 0.1 
and therefore N ^ 2/{—(3ia) ^ 4, which by 
no means is a large number; thus, one could 
question the large-A^ approximations used 
to obtain the result^ However, the impor- 
tant point is that one is in fact assuming 
that \K^a^\ > \Kia'^\ > |Kjva^|, i.e. the 
so-called ultraviolet-renormalon effect is by 
definition smaller than the usual perturba- 
tive error. It only gives a crude estimate 
of the maximum accuracy which can be ob- 
tained by computing higher-orders in per- 
turbation theory. It is not a new source of 
uncertainty, but just a lower bound on the 
same perturbative error. 

Any possible ultraviolet renormalon ambiguity 
is therefore already included in the perturbative 
error quoted in Eq. (38). Taking an extremely 
conservative attitude, one could argue that the 
onset of the asymptotic behaviour of the pertur- 
bative series has been perhaps reached already 
at = 3. Although, there is no signal of such 
behaviour in the presently known terms, the pos- 
sibility of a large TV = 4 contribution cannot be 
excluded. If that were the case, one should then 
take the size of the known TV = 3 term (and not 
a rough large- TV approximation to it!) as the ul- 
timate uncertainty of the truncated perturbative 
calculation. As explicitly shown in Eq. (|39|), the 
conservative perturbative error (|3^) we have been 
adopting [|| in the actual Rr calculation, already 
includes such a pathological possibility. 

4.3. Charm corrections 

Since the r mass is below the charm-production 
threshold, the calculation of Rr has been per- 
formed in the effective QCD theory with only 3 
active quark flavours. The decoupling of heavy 
quarks from the non-singlet currents p9| ] guaran- 
tees that the leading (logarithmic) charm-quark 
effects can be reabsorbed into the QCD run- 
ning coupling Jsot . These contributions are then 
taken into account, through the matching rela- 



® Moreover, a recent estimate [ [l8[ of the contributions 
coming from a double renormalon chain (i.e. two dressed 
gluon propagators) has shown that these higher-order 
effects are not suppressed, which makes the standard 
ultraviolet-renormalon calculus ill-defined. 
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tion between the values of as in the effective the- 
ories with 3 and 4 flavours |^.^. Nevertheless, 
there are additional corrections suppressed by in- 
verse powers of the heavy-quark mass, which a 
priori could be sizeable. These effects appear 
first at 0{a^), through a virtual c-c vacuum- 
polarization contribution to the gluon propaga- 
tor. The charm-quark contribution is fully known 
at 0(a2) and the O [a^ (M^/mc)"] con- 

tributions with n — 1,2 and 3 have been recently 
computed J53] . The final numerical effect is very 
small |52-|54|1: 



V TT / mi 
0.0004 . 



1 , f< 
225 ^°Hl^ 



107 

4500 



(45) 



4.4. Non-perturbative corrections 

The experimental data has already shown that 
the total non-perturbative correction to Rr is 
small [see Eq. (^6|)], in agreement with the theo- 
retical prediction. Nevertheless, with the present 
experimental accuracy, inverse power corrections 
introduce an uncertainty of about 1%. 

On the theoretical side, it is clear that those 
corrections are very suppressed (either by a fac- 
tor or a^/M^); however, the actual esti- 
mate of the leading contribution in Eq. (|2|) re- 
lies on phenomenological fits. The main source 
of information has been up to now e~^e~ data, 
which only measures ^.^Jy. In order to esti- 
mate the axial contribution, one needs to as- 
sume S'lj^/S^j^y ~ —11/7, as predicted by fac- 
torization. While the opposite sign of the vec- 
tor and axial-vector D — 6 contributions (helic- 
ity flip) does not depend on that approximation, 
the magnitude of their ratio is changed by non- 
factorizable corrections. Those non-factorizable 
radiative contributions have been found to be 
small [|5| (in the AIS scheme with a naive anti- 
commuting 75). Nevertheless, it would be desir- 
able to have a direct experimental measurement. 

The preliminary fits of the r-decay hadronic 
distribution, given in Table I result in a value of 
(^(6) which has the predicted magnitude but the 
wrong sign. However, given the present size of the 
experimental error bars and the strong correlation 
between the different fit parameters, one can only 



conclude that J'-^-' is in fact small. 

From the difference between the measured 
semi-inclusive r hadronic widths in Eq. (|3l|), 
Rr.y — Rt,a = 0.10±0.07, it is possible to extract 
the value of <5yl^. Note that the perturbative 
contributions cancel in this difference. Taking 
into account the small D = 2 and D = A correc- 
tions [Q, one gets ^yl^ — (6±5)%, in good agree- 
ment with the theoretical estimate in Eq. (|2l|). 
However, the present errors are too large for this 
determination to be significant. Clearly, better 
data is required. 

There is another source of non-perturbative 
contributions, which has been recently analyzed: 
instantons [^,g. It is well-known that instan- 
tons do spoil the OPE by introducing corrections 
that are power-suppressed by 9 or more inverse 
powers of the momentum. Owing to this large 
power suppression, the instanton contributions 
are very sensitive to the momentum scale: they 
are either negligible, and therefore irrelevant for 
practical applications, or very large, in which case 
they destroy the convergence of the OPE. Thus, 
instantons just provide a lower limit to the mo- 
mentum scale at which the OPE can be applied. 

The existing analyses |5^,^^ indicate that the 
instanton contribution to Rr is tiny. Unfortu- 
nately, it is difficult to make reliable quantitative 
predictions. For as{M^) = 0.32, the instanton 



effect could reach the 0.2-0.3% level fSTf, but the 
theoretical estimates are no longer valid at larger 
as values. Nevertheless, the instanton contribu- 
tions to Rt.v and Rt,a are predicted to be much 
larger than the total correction to Rr |0 (there is 
a cancellation between the vector and axial- vector 
contributions to Rr, similar to the one occurring 
m (5(6)). The possible size of the instanton effect 
can then be bounded using the present informa- 



tion on Rr 



y/A- 



The data points in Figure ^ have a very soft 
dependence on the energy scale, which discards 
any sizable instanton contribution above 1.2 GeV. 
Moreover, a recent phenomenological analysis of 
i?r,v — Rt,a, where the instanton effect is sup- 
posed to be maximal, concludes that the in- 
stanton contribution to Rr is smaller then 0.0005. 
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4.5. Uncertainties in the running of 

The final determination of as{M^) has to be 
compared with the values of as measured at other 
scales, for instance as(M|). The correct way of 
running the strong coupling through the differ- 
ent physical thresholds (c, b, . . . ) at C(a^) was 
known already twelve years ago [ pl[ . One con- 
siders two different strong couplings, al and 
as , in the effective QCD theories with Nf 
and Nf-i active flavours, and imposes a series 
of matching conditions [^l|j5^ ] which guarantee 
that the same physical predictions are obtained 
in both theories. Contrary to common belief, the 
MS strong coupling is (in general) not continuous 
at the matching point . The exact point where 
the matching between both theories is performed 
is actually irrelevant [ p9[ , provided one does not 
choose a crazy ii/mq(jjr) ratio which would in- 
duce huge higher-order corrections. 

Taking as{M'^) = 0.33 as input, the resulting 
value of as(M|) changes by less than 0.5% if the 
^-threshold matching point is varied between 2 
and 40 GeV, or if one changes the c-threshold 
between Mr and mi,. Note that the usual incor- 
rect procedure p5| of requiring a continuous as 
largely overestimates this error. 

The matching conditions take correctly into ac- 
count the logarithmic heavy-quark corrections. 
The physical observables get in addition l/m^ 
contributions that need to be computed. It has 
been advocated 1 45 4^ that the uncertainty asso- 



plicit calculation of each physical observable, as 
discussed in Sect. 4.3, The induced correction to 



ciated with these quark-mass effects can be esti- 
mated through a generous variation of the match- 
ing point. This is not quite correct. Varying 
the matching point, one is only testing the /i- 
dependence of the logarithmic corrections. As- 
suming that one knows the matching conditions 
to a sufficiently high order, the variation of /i 
would not produce any numerical change (i.e. 
zero error!); but heavy-quark effects would be of 
course still there. 

In fact, within the MS scheme, the l/rn^ 
contributions have nothing to do with the run- 
ning of as ■ The renormalization group equations 
only involve logarithms, because MS is a mass- 
independent renormalization scheme. Mass ef- 
fects are of course taken into account in the ex- 



Rt, given in Eq. ([451), is very small. 

5. SUMMARY 

Because of its inclusive nature, the total 
hadronic width of the r can be rigorously com- 
puted within QCD. One only needs to study 
two-point correlation functions for the vector and 
axial- vector currents. As shown in Eq. (^, this 
information is only needed in the complex plane, 
away from the time-like axis; the dangerous re- 
gion near the physical cut does not contribute at 
all to the result, because of the phase-space factor 
(1 — s/M'^Y- The uncertainties of the theoreti- 
cal predictions are then quite small. Notice that 
the accuracy of the Rr predictions is much better 
than the corresponding estimate of i?e+e-('S) at 
s = M^ [Re+e-{s) measures the vector spectral 
function on the physical cut, where the theoreti- 
cal predictions (at least at low energies) are more 
uncertain] . 

The ratio Rr is very sensitive to the value of 
the strong coupling, and therefore can be used to 
measure as{M^) This observation has trig- 
gered an ongoing effort to improve the knowl- 
edge of i?r from both the experimental and the 
theoretical sides. The fact that Mr is a quite 
low energy-scale (i.e. that as(M^) is big), but 
still large enough to allow a perturbative analy- 
sis, makes Rr an ideal observable to determine 
the QCD coupling. Moreover, since the error of 
as(/i^) shrinks as /i increases, the good accuracy 
of the Rr determination of as (M^) implies a very 
precise value of as{M'^). 

The theoretical analysis of Rr has reached a 
very mature level. Many different sources of pos- 
sible perturbative and non-perturbative contribu- 
tions have been analyzed in detail. As shown in 
the previous section, the final theoretical uncer- 
tainty is small and has been adequately taken 
into account in the final as{M'^) determination 
in Eq. (29). (Note, however, that there are still 
sizeable fluctuations in the Bi measurements , 
which could slightly modify the final result). 

The comparison of the theoretical predictions 
with the experimental data shows a successful and 
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consistent picture. The as(M^) determination is 
in excellent agreement with the measurements at 
the Z-mass scale, providing clear evidence of the 
running of as ■ Moreover, the analysis of the semi- 
inclusive components of the r hadronic width, 
Rt.v, Rt.a and i?r,s, and the invariant-mass dis- 
tribution of the final decay products gives a nice 
confirmation of the reliability of the theoretical 
framework. 
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